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^ ! Abstract 

The mergings of energy levels associated with the breaking of VT symmetry in the model 
of Bender and Boettcher, and in its generalisation to incorporate a centrifugal term, are 
analysed in detail. Even though conventional WKB techniques fail, it is shown how 
the ODE/IM correspondence can be used to obtain a systematic approximation scheme 
which captures all previously-observed features. Nonperturbative effects turn out to play 
a crucial role, governing the behaviour of almost all levels once the symmetry-breaking 
transition has been passed. In addition, a novel treatment of the radial Schrodinger 
equation is used to recover the values of local and non-local conserved charges in the 
related integrable quantum field theories, without any need for resummation even when 
the angular momentum is nonzero. 



1 Introduction 



Following a paper by Bender and Boettcher pP, itself inspired by a conjecture of Bessis 
and Zinn- Justin [2], the subject of "PT-symmetric quantum mechanics has attracted growing 
interest. (A couple of earlier papers exploring related themes are [31 E], while 0-112] provide 
a small sample of subsequent work.) Many studies have discussed general interpretational 
questions, but there are also mathematical issues to be addressed. One such is the topic of this 
paper: the pattern of "PT-symmetry breaking in the original model of Bender and Boettcher 
^ [3] , and in the generalisation of their model that was introduced in [5] . Numerical work has 
demonstrated an intriguing pattern of merging energy levels associated with the symmetry- 
breaking transition, but a complete analytic understanding has proved elusive, partly because 
standard WKB techniques break down in the regime where the symmetry-breaking occurs [Q. 
In this paper we show how this problem can be overcome, using a recently-discovered link 
with the theory of integrable quantum field theory [21] (the 'ODE-IM correspondence') to 
develop an approximation scheme for the energy levels which captures all features of the 
previously-observed behaviour in a controlled manner. 

The remainder of this paper is organised as follows. Section describes the symmetry- 
breaking transition in more detail, while section [3] discusses earlier theoretical treatments. 
The relevant connections with integrable quantum field theory are recalled in section 0J and 
then applied to the problem at hand in section Sectional then shows how a novel treatment 
of the radial anharmonic oscillator can be used to recover some of the results used earlier, 
calculations which may be of independent interest for the ODE-IM correspondence since they 
relate to the values of conserved charges in certain integrable quantum field theories. Finally, 
section contains our conclusions. 

2 Some VT phenomenology 

In the early 1990s, Bessis and Zinn-Justin 2 conjectured that the non-Hermitian Hamiltonian 

H = p 2 + ix 3 (2.1) 

should have a real and positive spectrum, if the boundary condition ifi £ L 2 (R) is imposed 
on the position-space wavefunctions ip{x). This initially-surprising proposal was motivated 
by considerations of the quantum field theory of the Lee- Yang model; it was subsequently 
put into a more directly quantum- mechanical context by Bender and Boettcher pQ, who also 
suggested that the spectral properties of the one-parameter family of Hamiltonians 

H M =p 2 -(ix) 2M (2.2) 

might be of interest, with M a positive real number. Varying M from 3/2 to 1 interpolates 
between H.3/2 , the Bessis-Zinn-Justin Hamiltonian, and H\ , the much more familiar Hamil- 
tonian of the simple harmonic oscillator. For M > 2, analytic continuation (in M) of the 
eigenvalue problem requires boundary conditions to be imposed on a suitably-chosen contour 
in the complex plane, away from the real axis pQ; given this, Bender and Boettcher found 
that the conjecture of Bessis and Zinn-Justin could be strengthened to the statement that the 
spectrum of TC M is real and positive for all M > 1. This reality property is associated with, 
but not completely explained by, the so-called VT symmetry of the corresponding spectral 
problems. In fact, it is only recently that a complete proof of the Bender-Boettcher conjecture 
has been given |10j . 
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Figure^] Real eigenvalues of p 2 — (ix) 2M + 1(1+1) /x 2 as functions of M, for various values of I. The 
(curved) dotted lines show the asymptotic locations of the level-mergings, to be discussed in sectional 
below. 

Most interesting, though, is the behaviour of the spectrum as M falls below 1, shown 
in figure Numerical results indicate that infinitely-many eigenvalues collide in pairs and 



2 



become complex, leaving a finite number of real eigenvalues [Q. As M decreases further, these 
successively pair off and become complex until finally the second and third levels merge and 
only the ground state remains real, which itself diverges to infinity as M — > (For M = \ , 
the spectrum of (|2.2j) is null, as can be seen by shifting ix to ix — E and solving the differential 
equation using the Airy function.) 

The mysterious nature of the transition at M = 1 is highlighted by a further generalisation 
of Bessis and Zinn-Justin's Hamiltonian, to incorporate an additional 'angular momentum' 
term |S]: 

n M , l = P 2 -(ixf M + 1 ^. (2.3) 

Here I is a real parameter, which can be assumed no smaller than —1/2, since the problem is 
unchanged by the replacement of I by — 1 — I . For I ^ 0, the contour on which the wavefunction 
is defined should be distorted below the singularity at x = , in addition to any distortions 
away from the real axis required for M > 2. The spectrum of TCm,i is then real for all M > 1, 
and positive if in addition I < M/2 |10j . As before, there is a transition at M = 1, with 
infinitely-many levels becoming complex. However, as illustrated in figure^, for sufficiently- 
negative values of I there is a remarkable change in the way that the remaining real levels 
pair off: the parity is reversed, so that the second level is paired not with the third level, 
but rather with the ground state, and so on up the spectrum. This reversed pairing holds 
true for sufficiently-high levels for all negative I, but as I gets closer to zero the level with 
which the ground state is paired moves up through the spectrum, leaving the connectivity of 
the original Bender-Boettcher problem in its wake and allowing for a continuous transition 
to the previous behaviour at I = 0. This is described at greater length in [S], and is best 
understood by looking at figures^, -Qi, or at figure 2 of [Hj; the challenge is to obtain an 
analytic understanding of why it occurs. 

3 Some VT Theory 

A first insight into the phenomena described in the last section comes from the observation 
that, while Hm is not Hermitian in any simple sense, it is invariant under the combined action 
of the operators V and T [Q, where V is parity and T time reversal, acting on Schrodinger 
potentials V{x) as 

VV(x)V- x = V{-x*) , TVix)?- 1 = V{x)* . (3.1) 

(The complex conjugation in the definition of parity ensures that the deformed contours 
required for M > 2 are mapped onto themselves, but is otherwise unimportant.) As shown 
in jS], VT invariance implies that eigenvalues are either real, or come in complex-conjugate 
pairs, much like the roots of a real polynomial (see also 20 ). Real eigenvalues correspond to 
wavefunctions symmetrical under VT, complex eigenvalues to a spontaneous breaking of this 
symmetry. In contrast to Hermiticty, on its own VT symmetry is not enough to prove reality; 
but it does mean that if a level is to become complex, then it must pair off with some other 
real level first. 

To deal with the transition at M = 1, Bender, Boettcher and Meisinger used a basis 
of the exactly-known eigenfunctions |n) of 7i\, the simple Harmonic oscillator: 

Hx\n) = (2n+l) \n) , n = 0,l,... (3.2) 
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and approximated the Hamiltonian for small e := 2(M— 1) as 



n 



l+e/2 



p 2 + x 2 + ex 2 [in \x\ + ^sgn(x)] + (9(e 2 ) . 



(3.3) 



Truncating the approximate Hamiltonian to the two-dimensional subspace spanned by |2n— 1) 
and |2n), for large n and small e the matrix elements are [H] 



H 



trunc 



o(2n-l) i6(2n) 
»6(2n) a(2n) 



where 



a(n) 



6 77* 

2n + 1 + — m(n) 



b(n) = - en 
V ; 3 



(3.4) 



(3.5) 



Diagonalising Htrunc yields the prediction that the energy levels of Tti+ e /2 should join and 
become complex for 

3 



8n 



(3.6) 



matching the numerical observation that the points at which levels merge approach the line 
M = 1 as the level n tends to infinity. 

However, the agreement of this result with the numerical data is only qualitative, and even 
at this level, there are problems. Firstly, the approximation also predicts a merging of levels 
for positive values of e, which is clearly wrong; and secondly, if one instead truncates to the 
subspace spanned by the levels |2n) and |2n+l), then the same approximation predicts that 
these two levels should pinch off. As underlined by the effect of the angular-momentum term, 
it is by no means obvious a priori how the pairing of levels should go, and it is reasonable to 
demand that any full understanding of the transition at M = 1 should include a robust (and 
correct!) prediction on this point. Figures |2Jt and 03 illustrate the situation. 
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Figure [2^: Truncation to levels 8 and Figure Eb : Truncation to levels 9 and 

9, compared with numerical data 10, compared with numerical data 

The discussion of j^j concerned the zero angular-momentum case. Since this is on the 
boundary between the two sorts of level connectivity, one might hope for a clearer signal 
from the truncated Hamiltonian at I ^ 0, where the matrix elements can be given in terms of 
hypergeometric functions. However, we have found that the problems described above persist. 
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The approximation (both for zero and for non-zero angular momentum) does improve if the 
dimension of the truncated subspace is increased, but this rapidly reverts to a numerical 
treatment, and any analytical understanding is lost. Alternatively, it might be hoped that 
the full set of matrix elements in the harmonic oscillator basis would simplify sufficiently in 
a suitable large-n, small-e limit that a precise prediction could be made. While this seems 
a technically-interesting avenue to explore, it is a delicate exercise to extract the relevant 
asymptotics, and we have yet to make useful progress in this direction. 

In more standard quantum-mechanical problems, methods such as truncation are not 
necessary when discussing the asymptotic behaviour of the high-lying levels, as the WKB 
approximation can be used. Indeed, this method gives very accurate results for the current 
problem when M > 1 pQ. However, the direct application of WKB fails for M < 1 T]: the 
path along which the phase-integral quantization condition should be taken crosses a cut and 
ceases to join the relevant two turning points. A naive analytic continuation of the WKB 
results for M > 1 to M < 1 also fails, producing energy levels which do not even become 
complex. 

In the following we shall show how these problems can be avoided, by developing an 
approximation scheme using the so-called ODE/IM correspondence, which links the Bender- 
Boettcher problem and its generalisations to a set of integrable quantum field theories asso- 
ciated with the sine-Gordon model. For M > 1 the WKB approximation is recovered, but 
for M < 1 the results are significantly different. In preparation, the next section gives a brief 
review of the correspondence. 



4 More VT Theory: the ODE/IM correspondence 

In its simplest form, the ODE/IM correspondence links the spectral properties of certain 
Schrodinger equations to functions which appear in the study of integrable quantum field 
theories in 1+1 dimensions. The phenomenon was first observed in |2l], for Schroodinger 
problems with homogeneous potentials V(x) = x 2M . The correspondence was generalised to 
incorporate an additional angular-momentum like term in and a more general potential 
in |26| , while its relevance to the non-Hermitian spectral problems of "PT-symmetric quantum 
mechanics was found in jSJ. A longer review of the subject is given in [27], and here we shall 
just summarise the key results, largely following jSj. 

Let {Ei} be the set of eigenvalues of TL M l , defined in (|2.3jl . and let {ej} be the eigenvalues 
of a similar-looking but Hermitian problem 

(-^ + >™)« = E, («, 

with so-called 'radial' boundary conditions ip{x) = 0(x l+1 ) as x — > , ip(x) — > as x — > oo. 
(Correspondingly, the earlier, PT-symmetric problems are sometimes said to have 'lateral' 
boundary conditions.) A pair of spectral determinants for these two problems, T{E) and 
Q(E), can be defined as follows. For any E, not necessarily an eigenvalue of either problem, 
let y(x, E, I) be the unique solution to (|4,lj) which has the same asymptotic as x —* +oo on 
the real axis as 

y WKB (x,E) :=^=P( X )-^ exp (- jf \/W) A (4.2) 
where P(x) = x 2M — E. (Note, so long as M is positive, I is not involved in the specification 
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of this asymptotic.) To give some examples: 



M>1: y(x,E,l) ~ _ x -W2 exp (__i_ x M + f) (43) 
M = l: y(x,£,Z) ~ exp(-±x 2 ) (4.4) 

l>M>i: y (x,iM) ~ ^i-^exp(-jfei M+1 + j4^) ( 4 - 5 ) 

More generally, the explicit asymptotic of y(x,E,l) changes whenever M = l/(2m— 1) with 
m a positive integer, as is easily seen from (|4.2[) . We also define tfj(x,E,l) to be the solution 
to Q4.1JI which behaves as x — > as 

i/>(x,E,l) ~x m + 0(x i+3 ). (4.6) 

Here we continue to take file I > —1/2, making the definition unique; if desired, a second 
solution ip(x, E, —1 — I) to 1)4.1)1 can be uniquely specified by analytic continuation in I. 

Denoting the Wronskian fg' — f'g of two functions f{x) and g{x) by W[f, g] , we now 
define T and Q by 

T(E) := W[y(iox,io- 2 E,l),y(u- 1 x,uj 2 E,l)] (4.7) 
Q(E) := W[y{x,E,l)^{x,E,l)} (4.8) 

where 

W = e W(M+l) _ (4>9) 

Then it is straightforwardly shown |5j that the zeroes of T are exactly the points {— E{\, and 
the zeroes of Q are the points {ej} , so these two functions are indeed spectral determinants. 
Furthermore they obey the following functional relation: 

T(E) Q(E) = uj- 1 - 1 ' 2 Q{oj- 2 E) + J +1 ' 2 Q{uj 2 E) (4.10) 

which is the same as Baxter's T-Q relation |28j from the theory of integrable lattice models, in 
the form which arises in connection with the quantum field theory of the massless sine-Gordon 
model (Strictly speaking, for I ^ —1/2, Q as defined here corresponds to the function 
denoted A in [22].) At the points M = l/(2m— 1) where the asymptotic form of y changes, a 
correction term appears in (|4.1U|) |291 |H] , but the only case potentially relevant in the following 
is M = m = 1, to which we shall return briefly below. 
Now define 

n (p\ . .21+1 Q(u 2 E) 

a(jB) - W Q(^E)- (411) 

By 1)4.10(1 . a(E) = —1 precisely at the zeroes of T(E) and Q(E), that is at the set of points 
{— Ei} U {ej} *. Furthermore, provided the {ej} lie on the positive real axis and the {— Ei} lie 
away from it, a(E) can be found by solving a nonlinear integral equation, as follows [3011311125] , 
Trade E for a 'rapidity' 9(E), such that 

e (m+i)/(2m) = re e (412) 



*Note, exceptionally, the zeroes of the LHS of l|4.10|l might coincide with simultaneous zeroes of the factors 
on the RHS, but this will not concern us below. 
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Then, for |3?m0| < min(7r, tt/M) , the 'counting function' 

f{6) :=\oga{E{9)) 

solves 



(4.13) 



f(Q) = iTr(l + 1) - imre e + / y(0-0') log(l + e s{8 ' ] ) dQ' 

-/ </>(0-0') log(l + e- / ^)d0 / (4.14) 



where 



and 



m 



7T ■ 



2M- 



M' 



-B 



3 1 ' 
2 ' 2M < 



r(p)r( 9 ) 



C t p - i (\ 

Jo 



t) q - l dt 



(4.15) 



(4.16) 



r(p + <z) 

is Euler's beta function. The integration contours C\ and C2 in (j4.14j) run from —00 to +00 
just below and just above the real axis, close enough that all of the points {0(ej)}, and none 
of the points {0{— Ei)}, lie between C\ and C2 ,, and the kernel function f(0) is 



<p(0) 



Jk9 



sinh(f^fc) dk 
2 cosh (p) sinh^fc) 2^ 



(4.17) 



The value of the normalisation factor r is arbitrary, but to match the conventions of |29l |H] 
one should take 



(2M+2)T 



M 
M+l 



(M+l)/M 



(4.18) 



The nonlinear integral equation ()4.14|) must be modified if any 9(ej) moves outside the 
contours C\ and C2 t 29 , or if any 9{—Ei) moves inside them [22j. However these possibilities 
do not arise in the current context, at least for small values of /: for M > 1 and |2Z+1| < M+l 
it can be shown that all of the ej lie on the positive real axis of the complex E plane, and all 
of the —Ei on the negative real axis (see (HI El)- F° r M < 1, the arguments showing that the 
{ej} lie on the positive real axis continue to hold, but, as we already saw, the {—Ej} do move 
away from the negative real axis and become complex. However, our numerical results show 
clearly that they never become near enough to the real axis to upset the arguments below. 



5 Asymptotics from the integral equation 

Equation (|4,14|) provides an effective way to solve the radial spectral problem (j4.1|) 24^ : given 
that the zeroes of T all lie away from the positive real E axis, any point on the real 6 axis at 
which f(9) = (2n+l)7ri for some n S Z must correspond to a zero of Q. The values of f(8) 
on C\ and C2 can be obtained numerically from (|4.14|) by a simple iterative procedure, after 
which the same equation can be used to obtain f{6) on the real axis, allowing those points 
at which f(9) = (2n+l)iri to be located with high accuracy. The leading approximation 
for f(9), found by dropping the integrals from (|4.14j) . yields the usual WKB result; higher 
corrections can be obtained using the asymptotic expansion for Q(E) in terms of local and 
non-local charges given in |29[ I33j . 
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The same numerical approach was used in 8 to analyse the non-Hermitian spectral prob- 
lems based on (|!2.3|) : it was also used to generate figure ^ above. However an analysis of the 
implications for the asymptotic behaviour of the energy levels was not given. Compared to 
the Hermitian case, there are a number of subtleties, which we now discuss. 

The first point, which was already taken into account in the numerical work described 
in [H], is that the zeroes of T(E) are either on or near to the negative real E-axis, which is 
the line Qm9 = vr(M+l)/(2M) on the complex fl-plane. Since (M+1)/(2M) > min(l,l/M) 
for all M 7^ 1, this means that the relevant values of f(9) lie outside the strip described by 
(|4.14|) . Instead, the so-called 'second determination' must be used [31] • This arises because 
the kernel tp{9) has, amongst others, poles at 9 = ±in and ±in/M, which add residue terms 
to the analytic continuation (|4.14|) outside the strip |9m#| < min(-7r, tt/M) . (See for 
related discussions in the context of the thermodynamic Bethe ansatz.) 

For M > 1 (sometimes called the attractive, or semiclassical, regime in the quantum 
field theory context) the poles in (p{9) nearest to the real axis are at 9 = ±iir/M, and have 
residue — Increasing ^sm9 past tt/M , the pole at 9 = —m/M crosses C\ and then C2, 
generating a term —f(9—iir/M) which should be added to the right-hand side of (|4. 14f) in 
order to give a correct representation of f(9) for 9 just greater than tt/M. Using the 
original representation ()4.14j) to rewrite the extra term, the resulting expression is 



f(0) = -i(l-e- in/M )mre e + / ip n (9-9 f ) log(l + e /(e,) ) dSf 

f ^ II (9-9')log(l + e- f(e ' ) )d9' (5.1) 

JCo 



where 



/7T 



2i cos^) sinh(0 - m 



VII (9) = m - ^(9-zn/M) = J™'^ V % • (5-2) 

7r(cosh(2#-f ) -cos(f )) 

An examination of the poles in ipjj(9) shows that (|5.1jl holds for tt/M < Qm9 < ir, which 
includes the neighbourhood of the line 9m 9 = n(M+l)/(2M) relevant to the hunt for the 
zeroes of T(E). (For simplicity we have supposed that C\ and C2 are infinitesimally close to 
the real axis; otherwise an 'intermediate determination' is also required, applying when the 
pole in (p{9) has crossed C\ but not C2 .) 

For M < 1 (the 'repulsive regime' in quantum field theory language) the nearest poles are 
instead at 9 = ±i7r, and for Qm9 > ir the extra term to be added to the right-hand side of 
(|4.14|) is equal to +f{0—m). Rewriting using the original representation as before, one finds 

f(9) = 2m{l+l)+ [ ip II (9-9')log(l + ef^)d9' 

JCi 

-[ fn(9-9')log(l + e- f ^)d9' (5.3) 
Jc 2 

where 

Msin(7rM) 

VI m = m + = <cosH2M6 _ mM) _ cos(7rM)) • (5-4) 

This holds for vr < 9m 9 < tt/M. 

The goal is to locate the zeroes of T(E) on or near the negative-S axis. To this end we 
set 9 = iTr(M-\-l)/ (2M) + 7, define the shifted counting function 

5(7) := f{m{M+l)/(2M) + 7) (5.5) 



S 



Abusing the notation a little, we shall sometimes write g as a function of the energy E, 
g{E) := f(-E), where 

E = ( re 7)2M/(M+l) (5 6) 

The eigenvalues of the 'PT-symmetric problem Q2.3|) are then those values of E such that 

g(E) = (2fc+l)vri (5.7) 

for some k £ Z . 

Defining ^(7) := (pu(iir(M+l)/(2M) + 7), 5(7) is given exactly, for |3rn7| < | l^jpH , 
by the following expressions: 

M > 1: 

5(7) = 2isin( 1 f ? )mre 7 + f ip(j-9') log(l + e /(9,) ) dfl' 

- / r/>( 7 -0 / )log(l + e- / ^'))^ / (5.8) 

with 

= 2cos(^)cosh(7) = f ( _ 1)fl+1 2 cos( ^ (2n _ 1)} e -( 2 n-l) 7 (5 _ g) 
n " 7r(cosh(2 7 )+cos(^)) ^ V > vr ^ " 1 ; 

M < 1: 

5(7) = 2ivr(/+i)+ / V(7-^)log(l + e /(9 ' ) )^' 

_ f ^(7-0 / )l O g(l + e - /(e ' ) )^ / (5.10) 

with 

/ n — Msin(-7rM) ,„2M , . 

^(7 = -7 w rr^ / = > (-1) sin TrMn) e - 2Mn T . (5.11) 

vyu vr(cosh(2M7) +cos(vrM)) ^-J V 7 vr v 7 v ; 

(We have recorded the large-7 expansions for ip as they will be needed later.) 

The formulae (|5.8[) and (|5.1()|) express 5(7) in terms of the values of f(0') on the contours C\ 
and C2, where the original ('first determination') specification ()4.14|) applies. For a given value 
of 7, the dominant contributions to the integrals in (|5.8j) and (|5.1()|) come when $te9' ~ 3ie7, 
since the kernel ^(7) is peaked about 7 = 0. For large ?R.e6' on C\ and C2, /(#') ~ — 
imre 9 . Since C\ is just below the real axis, and C2 just above, in this limit log(l + e^ e ') —* 
on Ci, and log(l + e~^ e >) —> on C2. Therefore, for large 3fte7 the integrals in (|5.8|) and 
(|5.10|) can be dropped, giving the following leading approximations^ for 5(7) : 



M > 1 
M = 1 
M < 1 



(M+l) 



3(7) ~2isin(2^)m£T- (5.12) 
5(7) =*7r(l+i) +ifS (5.13) 

9(l) ~ 2ivr(/+i) (5.14) 

^The result for M = 1, which is exact, follows directly from 14.141 since <p(6) = in this case and no second 
determination is required. Since m| M=1 = 7r/2, it interpolates between the results for M = 1 + and M = 1". 
This interpolation will be discussed in more detail around equation l|5.81|l below. 
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For M > 1, the leading approximation for the energy levels is obtained by setting the RHS 
of (fS"T2|) equal to (2n+l)7ri . Then, using (|4~TI)() . the WKB result of [TJEj is recovered: 



E n 



'3 + -L- 



2M/(M+1) 



sin(^)r(l+^) 



(5.15) 



For M = 1, there are some subtleties. Since the eigenvalues at this point are easily found 
exactly, they can be ignored, but their resolution is perhaps interesting. At M = 1, u = i 
and the factors Q(lo ±2 ) in Q4.11JI cancel, making / and g as naively defined by 1)4.13(1 and 
1)5.5(1 equal to n(l + \). However this is too quick: at M = 1 the TQ relation 1)4.10(1 is 
'renormalised' to T(E) Q(E) = ^-i/2+£/2 Q( W -2_ E ) + J+1/2-E/2 q^E) H3IB], and as a 
result ()4.11(l is naturally replaced by a(E) = uj 21+1 ~ e Q(ui 2 E)/Q(uj~ 2 E), which indeed matches 
(UTTTfl) . Next, the renormalised TQ relation implies T{E)Q(E) = 2cos({2l+l-E)n /4)Q{-E). 
Knowing that the zeroes of Q are positive and those of T negative allows this relation to be 
disentangled: the zeroes of Q must be at E = 21 — 1 + 4n, and those of T at E = 21 + 3 — An 
and —21 + 1 — 4n, for n = 1,2, ... . Recalling that the eigenvalues of the lateral problem 
are at the negated zeroes of T recovers the spectrum of the 'PT-symmetric simple harmonic 
oscillator', previously obtained in |H] and jS] via explicit solutions of the differential equation. 

Finally, for M < 1 the leading approximation to 5(7) is a constant, and does not give 
any information about the energy levels at all. This gives a novel insight into the observation 
of that standard WKB techniques fail for M < 1 : from the point of view of the ODE/IM 
correspondence and its associated nonlinear integral equation, it can be traced to the change 
in the nature of the second determination in the repulsive regime. 
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2 4 6 



10 12 14 




2 4 6 



12 14 



Et: M = 0.8 



Et>: M = 1 



Oi: M = 1.2 



FigureEl 3ro g(E) from the nonlinear integral equation. In all three cases, I = —0.001. 
The horizontal lines show the 'quantisation levels' g = (2n+l)7ri. 

To proceed further, the more detailed behaviour of g(E) is required. Figure shows 
g(E) in the various regimes, the plots being obtained by solving the NLIE numerically, 
but with no other approximations. Note that the nature of the energy level quantisation is 
radically changed for M < 1: even for the levels which are still real, there is no natural way 
to allot them a unique WKB (Bohr-Sommerfeld) quantum number via the counting function, 
as happens for standard eigenvalue problems in quantum mechanics, and also for these lateral 
('PT-symmetric) problems when M > 1. 
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The next task is to understand figure El analytically. It turns out that the direct use of the 
asymptotic expansion for Q(E) developed in OH] fails to capture most of the important 
features of the energy levels for M < 1 : roughly put, the wiggles in figure OK ar e missed, and 
the asymptotic behaviour of just a single level is captured - the partner of the ground state, 
which moves up to infinity as I — > 0~ (this can be seen on figure^, and is illustrated in more 
detail in figure below). The remaining levels only emerge when a 'beyond all orders' term 
is added. The need for this term is most easily understood if we continue to work directly 
with the nonlinear integral equation. 

Higher corrections to 5(7) can be obtained through a steepest-descent treatment of the 
integrals in or (|5.1t)j) . expanding the kernels ip^j — 6') for large 7 — 8' and using again 
the fact that the terms log(l + e^ e )) and log(l + e~^ 9 )) decay to zero as 5?e#' — > 00 along 
C\ and C2 respectively. This requires the contours C\ and C2 to be shifted so as to maximise 
the rates of decay of log(l + ') along them. For large Ke 6' with \Qm 6'\ < min(7r, ir/M) , 
f(O') ~ Z7r(/+|) —imre 6 , so C\ should be shifted down to Sm 6' = —ir/2 , and C2 shifted up to 
Qm6' = +7r/2 . (For M < 2, these lines do indeed lie inside the first determination, for which 
the just-mentioned asymptotic of f(9') holds, and we restrict to such cases in the following.) 
These shifts of contours are not entirely innocent operations, as the kernel functions 
have poles at 7 = ±^ ■ Taking C\ and C2 past these poles generates residue terms 
proportional to the values of log(l + e^^ 9 **) there. These are the 'beyond all orders' terms 
just advertised. 

Explicitly, let C\ and C2 be the contours C\ and C2 shifted onto the steepest descent (for 
large ReO') paths Qm6' = —tt/2 and ^smO' = +7r/2 respectively. For 1 < M < 2 , (|5,8|) can 
then be rewritten as 

5(7) = 2isin( 2^) mre 7 

+ [ V(7-^)log(l + e /(9,) )^'- [ V(7-^)log(l + e~ /(e,) )^' 
JCi Jc 2 

+ log(l + e /(7_i ?^ H) ) - log(l + e-^+TT 1 ^)) . (5.16) 

For g < M < 1, the representation (|5.10[) is rewritten in a similar manner, modulo a change 
in the signs of the extra terms, due to the opposite residues of the relevant poles. (For M < i, 
the poles are not encountered in the shifting of C\ and C2, and so no extra terms are needed.) 

Once the contours have been shifted, the kernel functions are expanded using (|5.9j) and 
(j5.11|) to give asymptotic expansions of the integrals in the standard way. Setting E = 
(re 7 ) 2M /( M+1 ), for 1< M < 2 this yields 

(A/+1) 

g{E) ~ 2ism(^ T )mE^M- 
00 „ , 

+ i J>ir +1 ^ cos(^(2n-l)) E-^^) - g nonpcvt (E) (5.17) 

?1=1 

where the coefficients b n are given in terms of f(8) by 

ib n = f {re e ) 2n ^log(l + e f(e) )de- f (re 9 ) 2 ^ 1 log(l + e~ /W ) d9 (5.18) 
JCi Jc 2 

and 

/l + e -/(7+f^)\ 

5nonpert(^) = log ~ | M -iK ■ ( 5 - 19 ) 
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(After the expansion has been made, the contours in l|5.18j) can be shifted back to their original 
locations without encountering any poles.) The geometrical interpretation of <? non p er t W1U be 
clearest if the branch of the logarithm in (|5,19|) is taken in the interval [—n,ir), and this will 
be assumed from now on; note that this choice has no effect on the quantization condition on 
the eigenvalues. 

For | < M < 1, the result is instead 

00 9 M 

g(E) ~ 2,vr(/+i) + i £(-l) n sin(vrMn) E'( M +^ + g DOQpcit (E) (5.20) 

n=l 

with coefficients c n given by 

f (re e ) 2Mn log(l + e^)de- [ (re e ) 2Mn log(l + e~M) d9 (5.21) 
JCi JCo 



%c r . 



and fl , n0 npert(-^) 1S again given by (|5.19f) . Note that <? nonpert now contributes with the opposite 
sign, while the modulus signs in the formula mean that the numerator is always evaluated 
above the real axis, and the denominator below, if 7 is real. 

The integrals for the coefficients b n and c n are known exactly in terms of the ground state 
eigenvalues of the local and non-local conserved charges of the zero-mass limit of the sine- 
Gordon model on a circle |361 129j . Translating* the results from these papers into the current 
normalisations, 



vr 3 / 2 (4M+4f r((a+sk)(2n-l)) 

i o i w 1 ^ 7v\ hn-l , [p. 11) 



" [ 2n ~ l r(^(2n-l)) 

cos(vrMn) H n . (5.23) 



2TT ° 2Mn 



M (M+l) 2 ™ 

Here, hn-i and H n are respectively the ground state eigenvalues of Ibn-i , the n th local 
conserved charge, and M n , the n th (dual) nonlocal conserved charge. The precise definitions 
of these charges can be found in [HEl EB E3] (see also [HZ1 EH] ) • Explicit expressions for the 
local charge eigenvalues can be found up to high order; the first three are 



h 



1 



(4M+4) 
1 

(4M+4) 2 
1 

(4M+4) 3 



A 



M+l 



6 

A 4 - (M+1)A 



o (M+l)(4M+3)(M-3) 
60 



5(M+1) 4 (M+l)(4M 2 -23M-23) 2 
~~2~ A 12 A 

(M+l)(96M 4 -340M 3 +85M 2 +850M+425) 
1512 



(5.24) 



where A = The nonlocal charge eigenvalues are harder to calculate, but the first one is 

also known in closed form: 

~ = (M+i) 2 r(M+i)r(-2M-i)r(M+i+A) 
1 r(-M)r(-M+A) ' 



*The dictionary [HEg is /3 2 = 1/(M+1), p = (2Z+l)/(4M+4) . 
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Finally, the beyond-all-orders term fl'nonpert 

(E) in (jfTPTj) and (|5~2T)|) should be treated, 
which contains the so-far unknown function f{9). We first note that for 7 real f{^— ^—jp^) = 
-/* (7 + f¥) ( the * denoting complex conjugation) and so 

Snon per t(£) = 2*arg (l + e -/(?+f . (5.26) 

Given the branch chosen for the logarithm just after (|5,19|) . the value of the argument in ()5.26|) 
lies between — ir and ir. In principle the full asymptotic expansion for /, to be given in the 
next section, could now be substituted into this formula. However the interesting structure is 
already seen if the leading behaviour is used, which, from (|4.14l) . is 

/(7 + f ^■))^m{l+\)-ime^ ]M ^E^. (5.27) 

Hence 

ff aonpert (^)~2iarg(l + pe^) (5.28) 

where 

p(E) = e-msm^^E^r. (5 _ 29) 

(f)(E) = -tt(Z+|) + mcos^^E^r . (5.30) 

The final formula can be given a pictorial interpretation: gnonpert(-E) is approximately equal 
to 2i\{E) , where x(E) the interior angle of the triangle shown in figure |IJ 




1 



Figure H The trigonometry of g n on P mt(E) . 

Notice that p{E) is always less than 1 for M/l, and that as E increases with M fixed 
the point 1 + pe 1 ^ executes a diminishing spiral about 1, giving x(E) a damped oscillation 
between —ir/2 and 7r/2 which captures the behaviour seen for M < 1 in figure |2K above. 
The nonperturbative nature of ffnonpert (E) is easily seen: |x(-^)| is bounded by ^p(E) ; and 
at fixed M, the expansion parameters for the perturbative series parts of (j5,17|) and ()5.20|) 
are e := E~^m~ and r] := E ,_ ( M+1 ) respectively, in terms of which ^p is equal to either 

7r / ■ \M—1\ / \ n I • IM-ll / 

2 exp(— man 1 M ' /£) or ^ exp(— msm M 1 /772MJ. 

If on the other hand M is taken very close to 1 with E remaining finite, the nonperturbative 
term plays a key role. In this limit, p{E) ~ 1 and hence, from figure |IJ x(E) ~ (p(E)/2, and 

Snon P crt(£) « + im E^ . (5.31) 

This is just the behaviour required if the nonperturbative term is to 'smooth' the discontinuous 
change in the leading asymptotic of g(E) as M passes through 1, seen in equations ()5.12|) 
- (|5.14|) above. The phenomenon is reminiscent of the smoothing of Stokes's discontinuities 
discussed in (30]. Strictly speaking, since x remains between —ir/2 and ir/2 the limit of 
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SWipert (-^) A as M -> 1 is a sawtooth function obtained by returning the imaginary part 
of the RHS of (|5.31|) . modulo 2tt, to the interval [— tt, it]. This has an important effect on 
the counting of eigenvalues. Consider for example the limit M — > 1~ at fixed E: since the 
coefficients in series part of (j5.2U|) are all of order e = 2 — 2M, the limiting form of g(E) is 



9(E) 



M=l~ 



27ri(i + 1; 



TTl 



[-1,1] 



(5.32) 



where the subscript on the last term indicates that it should be returned, modulo 2, to the 
interval [—1,1] - so [x][_ L1 ] := (x+l)mod2 — 1. Since g is continuous for all M < 1, its 
M — ► 1~ limit has segments of infinite gradient at the points I + ^ — ^E = (2/c+l) where 
the corresponding sawtooth function would have a discontinuity. The quantisation condition 
g(E) = (2n+l)iri is therefore always satisfied at these values of E, in addition to the points 
I + ^ + ^E = (2k+l) which would have been found even before taking the sawtooth effect into 
account. Taken together, these points match the full spectrum of the PT-symmetric simple 
harmonic oscillator, discussed from another point of view just after equation (|5. 15f) above; the 
crucial role played by the nonperturbative term in this process is interesting. 






.7 0.8 0.9 1 
M 



Perturbative term 



|SJd: Nonperturbative term 



Et: Both terms 



Figure [SJ Different approximation schemes for the energy levels (small crosses) compared 
with the exact levels (continuous lines), for I — —0.001. 

More generally, for M < 1 and / ^ the nonperturbative term explains the infinite se- 
quence of merging levels high in the spectrum, while perturbative term is responsible for the 
single level joined to the ground state which moves upwards as / — > 0. To understand the 
reversed connectivity below this level, both effects must be incorporated simultaneously. Fig- 
ure [5] illustrates this by comparing the exact energy levels, computed using the full nonlinear 
integral equation, with the levels obtained by combining the leading asymptotic of g(E), (|5.12|) 
- (j5.14|) . with either a) the first perturbative correction, or b) the leading approximation to 
the nonperturbative correction, or c) both. (Note that the lowest group of crosses in figure 
[S^ at intermediate values of M is spurious, being caused by the over-large contribution of the 
perturbative correction at such small values of E.) 
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Explicitly, the crosses on figure mark points at which gapprox(-E') = (2k+l)iri, where 

2nonpcrt(£) iM > 1); 



5approx (E) 



° 1 V 2Af ^ V 2 2M ' 



i7T 3 / 2 r(M+f+Z) 



4iwyi-r 2 J f r(A/+ | )r( _ M)r( _ M+ i + ^ & t- </nonpert 



(£0 



(M< 1), 
(5.33) 

with fl , nonpert(-E') given by the leading approximation (|5.28j) . and m by (|4.15|) . For figure 
only the first and second terms on each RHS were used for ^approx! for figured, only the first 
and third. 

Even for the low-lying levels, the combined approximation does very well. Figure H3 shows 
a more stringent test by zooming in on the first three level-mergings; again, the match is very 
good. The errors are greatest in figure HJ}, near to the second merging. This is not surprising 
since the recombination of levels occuring there makes the locations of the eigenvalues partic- 
ularly sensitive to errors in the value of g(E). These errors become smaller for recombinations 
higher up the spectrum. 



5.5 



3.5- 
3 

2.5 




0.66 0.7 0.74 0.78 0.82 

M 

Et: First merging 




).66 0.7 0.74 0.78 0.82 
M 

Eb : Second merging 




0.84 0.88 . 92 

M 

Et: Third merging 



Figure The perturbative-plus-nonperturbative approximation for I = —0.001 (small 
crosses) compared with the exact levels (continuous lines) near the first three level mergings 
shown on figure [SJ:. 

To analyse the level-mergings for M < 1 in more detail, we take e = 2M—2 small and 
negative, and E fixed but large. Then the limiting form of the quantisation condition arising 
from (|5.33|) can be written as 



7T 



x (E) = nk-7rl + -(l+±)(l 2 -±)\e\E- 



(5.34) 



where fceZ, x{E) = ^ ffnonpert (E) is the angle shown on figure 01 and, to the same approxi- 
mation, the other quantities on the figure are 



p(E) 
HE) 



e a 



-*{l+\) + \E + ^\e\EHE). 



(5.35) 
(5.36) 



Since x(E) is always between —tt/2 and 7r/2, all solutions to (|5.34|) for I small have k = 0. If/ 
is nonzero, then sufficiently high in the spectrum the final term on the RHS of (|5.34j) becomes 
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insignificant and the condition reduces to x(-E') = —ttI. This is illustrated in figured given 
that e is small, p behaves as a slow mode in E, and as a fast mode, and so the point 1 + p e 1 ^ 
moves around the circle of approximately-constant radius p, this radius slowly decreasing as 
E increases. Eigenvalues occur when this point crosses the ray from the origin with argument 
—irl. Two such rays are shown on the figure: one, the upper dashed line, is for I < 0, and the 
other, the lower dashed line, is for I > 0. 



L _ I < 
/>(> 



Figure The approximate quantisation condition for 1^0. 

If p becomes too small, the circle and the ray cease to intersect and so the eigenvalue 
condition can never be met, no matter what (real) value is taken by the fast mode <j) . This 
means that the eigenvalues have merged to become complex. The critical value of p at which 
this first occurs is | sin(7rZ)| ; hence, the merging of eigenvalues occurs asymptotically along 
the curve 

E = _ 8 H*°frOI • (5-37) 
it \e\ 

The curved dotted lines on figures [Ijp and ^1 allow Q5.37|) to be compared with the 
behaviour of the exact levels. 

We can also predict which eigenvalues pair off, by thinking about the fast mode (p. If E is 
kept approximately fixed and |e| is increased (since e < 0, this corresponds to decreasing M), 
then from (|5.35|) p decreases while <p is unchanged. This makes it easy to see that it is pairs 
of eigenvalues associated with neighbouring crossings of the ray arg(z) = — ttI on the same 
side of the real axis which merge. Now as E increases from zero, the phase increases from 
— 7r/2 — 7rZ, meaning that the point 1 + pe 1 ^ starts below the real axis. If I < 0, it first crosses 
into the upper half plane before encountering the ray arg(z) = — irl, which it then does twice 
before crossing the real axis again. Hence the lowest level is paired with the second level, the 
third with the fourth and so on. If on the other hand I > 0, then starting again from E = 
the point 1 + pe 1 ^ encounters the ray arg(z) = —ttI once, below the real axis, then moves 
anticlockwise through the upper half z-plane before the next two crossings of the ray, both 
below the real axis, occur. Hence the connectivity of the levels is reversed, with the second 
and third, fourth and fifth, and so on, levels pairing off. This is exactly as observed in our 
earlier numerical work, now predicted analytically. 

Strictly speaking, we have just found the connectivities that the levels would have were 
they quantised by the nonperturbative term alone, since we discarded the final term on the 
RHS of (|5.34[) before commencing the analysis. However, if E is taken large enough before e is 
sent to zero, then this term dominates and so our argument also predicts the connectivity of 
the true level-mergings, high up the spectrum. Lower down, the behaviour of the true levels 
may differ from that of the nonperturbative-only approximation, as can be seen by comparing 
figures |Hb and Et- This corresponds to the recombination of levels, and simple estimate of 
where it occurs can be obtained by finding the intersection of the curve (|5.37|) with the curve 
^(Z+|)(/ 2 — j) |e|-E -2 = 7rZ, which approximates the energy level found by quantizing according 
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to the perturbative term alone, shown for I = 



—0.001 in figureEH. For small I < 0, this yields 



... /In |ttZ| \ 1/3 
Erecomb(') — I ^ 2 ; J ' (5.38) 

For I = —0.001 , -Erecomb ~ 8.36, which matches well the behaviour shown in figures^ and[3 
It remains to treat the marginal case, I = 0. This is much more delicate. Easiest to 
discuss is the connectivity of the levels. The last term in 1)5. 34|) functions much as the term 
— irl did before from the point of view of the fast mode <f>; given its sign, this means that 
the connectivity of levels for I = is the same as that for I > 0, which is indeed what was 
observed numerically. (A similar argument explains why the connectivity is reversed below 
E = -EVecomb for I < 0.) To estimate where the level merging occurs, note that the critical 
value of p is now sin(^|e|E -2 ). For large E and small e we must therefore solve 

e-xH 2 * = §|e|£T 2 (5.39) 

to find the critical value of E. The solution can be expressed in terms of the Lambert W 
function W(x), which is defined to satisfy W(x)e w ( x ' = x. (More information on this function 
can be found in 0Tj; it is related to the tree function T(x) by W(x) = —T(—x).) Explicitly, 

16 / ^5/2, ,3/2 

E = —5TT W -A- J I ) ( 5 - 4 °) 



vr 2 |e| V 32\/2 

where the subscript —1 selects the relevant branch of the W function (see |U). Now using 
the asymptotic T-i(x) = — W-\{— x) ~ ln(x _1 ) + ln(ln(x -1 )) for x — > + , the curve 1)5.3711 . 
which diverges as I — > 0, should be replaced by 



E. 16 



7T 2 |e| 



In | 3 .yLWl„h/ 32V ~ 2 



7T 5 /2| e |3/2y Vvr 5 / 2 |e| 3 / 2 



(5.41) 



This corrects the result of 0, equation ()3.6)1 above, and shows that while the two-by-two 
truncation used in that paper gives a qualitative understanding of the situation, is too crude 
to capture precise asymptotics. However, even with the second logarithm included in (|5.41|) . or 
using the W function directly, the approach of the level-mergings to their asymptotic locations 
for I = is much slower than that found for nonzero values of I. The dotted line on figure^ 
compares ()5.4U|) with the exact levels, obtained by solving the nonlinear integral equation. 



6 Asymptotics for the radial problem 

The approach of the previous sections can also be applied to the more-standard radial problem 
(|4.1j) . The analysis is very similar to that already given, though simpler because the subtleties 
of the second determination do not apply. 

The large-0 expansion of the first-determination kernel 1)4. 17|) is 

oo oo 

<P(0) ~ c ot(jt(2n-l)) e - {2n - l)e + ^ — tan(^Mn) e~ 2Mnd (6.1) 

n=l n=l 
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and so for large real 6 , 



CO . 

bn , „ ,„ 



/(E) ~ ^(/ + i)-i^^ co t(^(2n-l)) J B- L w^(2n-i) 

n=0 

00 . , 

+ i J2 — tan(vrMn)£-( M+1 )™ , (6.2) 

n=l 

where b n , c n and m are as before, and we denned 60 = ~ 7r ^ an (jM) m so as to absorb the 
exponentially-growing part of the asymptotic into the first series. The condition for E to be 
an eigenvalue of the radial problem can be written as 

f(E) = -(2k+l)iri, k £ Z (6.3) 

Written directly in terms of the conserved charges, and setting i_i := 1, the expansion 
coefficients of / are 

1 ^2 ' l 2M J 

Mc 2 2Mn ~ 

i — tan(vrMn) = -2% sm(TrMn)— — ^ # n . (6.5) 

Notice that (|6.2|) is a double series, unlike the single series (|5.17j) and (|5.2U|) found for 
the lateral problems, and that its form does not change as we move between the attractive 
(M > 1) and the repulsive (M < 1) regimes. For M < 2 the steepest-descent contour 
lies inside the first determination, at least for large E, and so, unlike for the PT-symmetric 
problem discussed above, there is no need to add a non-perturbative term. The double series 
could also have been recovered directly from the asymptotics of Q(0) found in [55], given the 
ODE/IM correspondence. Taking suitable care about the branch choices, the same applies 
to the series parts of the expansions for the lateral problems obtained in the last section, but 
the beyond all orders contributions would have been missed. 

For the radial problem, the WKB method should work for all values of M, and it is natural 
to try to recover the asymptotics from a direct analysis of the differential equation. In the 
rest of this section we show how the contributions to (|6.2j) from local conserved quantities can 
be found using an all-orders WKB approximation, and give some hints as to the origins of the 
nonlocal parts. Even for the local part of the discussion, there are some interesting features, 
as the application of the WKB technique to radial problems is tricky |421 1431 144j . 

A WKB treatment of the radial problem might begin by writing (|4.1j) as 

'' 2 ~ \ = 0, Q(x) = x 2M - E + . (6.6) 



x' 2 



For / and E positive, Q(x) has two simple zeroes on the positive x axis: one, x$, near x = 
and the other, x\, near x = E l / 2M . The leading WKB quantisation condition would then be 

^ Ixo V ~Q( X ) dx = (2k+l)n , k G Z . However, this does not give good results. The reason, 
first stressed by Langer [15], is that the WKB approximation derived from (|6.6|) breaks down 
near x = 0. Instead, Langer suggested making a preliminary change of variable 

x = e z (6.7) 
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and gauge transformation 



so that (j6.6j) becomes 



y(x) = e z / 2 ip(z) 



$ = , Q(z) = e {2M+2 *> z - Ee 2z + (1+1/2) 



The quantisation condition for the modified problem is 



-Q0) dz = (2/c+l)vr , k£ 



(6.8) 



(6.9) 



(6.10) 



which, changing the variable back to Inz = x, amounts to the trading of 1(1+1) with A 2 = 
(I + 1/2) 2 in the original potential. This is sometimes called the Langer modification, and it 
allows the energy levels of the radial harmonic oscillator and the radial Coulomb potential to be 
recovered exactly. However, away from these points, the generalisation to incorporate higher- 
order WKB corrections is difficult. In particular, the leading Z-dependence in (|6.2|) . namely 
the term iTr(l+^), only emerges once an all-orders resummation has been performed |44j . 

The problem is that I and E are badly tangled up in formulae such as (|6.1U|) . To get 
round this, we take the Langer-transformed equation ()6.9j) and change variable z —* jz. 
Transforming back to the variable x, using (|6.7|) and (|6.8|) in reverse, we find 



dx 2 



+ 7 x 



2 f2M+2) 7 -2 



E, 2 x 2 ^ 2 + ^ 



y = 



(6.11) 



where 



Z = 7 A-l/2. (6.12) 
If 7 = 1/(2A) then [=0, and, changing variable to t = E~^' M x, (|6. llf) simplifies to 



dt 2 



+ Q(t) 



where 



Q(t) 



4A 2 



l/\-2, t M/X _ ^ 



y(t) = 



E 



-(M+1)/2M 



(6.13) 



(6.14) 



A key feature of this equation is that the E'-dependence, contained in e, has been factored out 
of the transformed potential Q(t). In fact, e takes the role of the variable which organises the 
whole higher-order WKB series, in eq. (|6.24j) below. This means that a systematic WKB treat- 
ment of the transformed equation should yield the asymptotic behaviour of the energy levels 
as E — > oo at fixed M and I, which is just what we want. The contrast with (|6.9|) accounts for 
the difficulties encountered when trying to work directly with the Langer-transformed equa- 
tion. The simplification has a price, though: for / 7^ the turning point which had been near 
the origin is replaced by a singularity of order 1/A — 2 exactly at x = 0. This requires special 
treatment, which we illustrate by calculating the leading correction explicitly by asymptotic 
matching. 

To build a solution which to leading ('physical optics') asymptotic approximation solves 
the transformed eigenvalue problem, the t axis can be split into four regions, defined for 
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e — > + as follows^: 



Region I 
Region II 
Region III 
Region IV 



t > 1, \t-l\ >e 2/3 
|t-l| < 1 

1 > t > 0, |t| > e 2A , > e 2/3 

i > 0, \t\ <C 1 



(6.15) 



In regions I and III, the physical optics approximation is good, while regions II and IV need 
separate treatment. In region I, the decaying WKB solution is 



1 



[Q(t)]V* GXP 



1 



\J Q(u) du 



(6.16) 



This can be continued through region II, which contains the simple turning point at t = 1, 
using an Airy function in the standard way [IS], to find y(t) ~ Vmit) in region III, where 



3/m(*) 



-Q(t)]V4 



cos 



1 Z" 1 7T 

-y 4 



(6.17) 



Taking t near zero, but still in region III, this behaves as 



2/m(*) 



■ cos 
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(6.18) 



This solution breaks down as t —* 0, and must be matched to a solution of the approximate 
ODE in region IV, which is 



d 2 



1 



dt 2 4X 2 s 2 

Imposing y — > as t — > 0, the appropriate solution is 

y w {t) = (3yftJ x am) 



y{t) = 



(6.19) 



(6.20) 



where J\ is a Bessel function, and j3 a so-far arbitrary normalisation constant. (The behaviour 
as jt2* — > is y IV (i) ~ /3t (2e)~ A /r(A+l), as follows from the £ — > asymptotic ~ 
(£/2) A /r(A+l).) The large-£ formula J A (f) ~ a/Jcos [£ - |tt - f ] then implies 



/x „ /2e ,1 1 / 1 1 A-7T 7T \ ,. \ 

Vivit) ~ P\ — t2 4A cos -t2A -] for i > e . 

V 7r V e 2 4 



(6.21) 



For E to be an eigenvalue, this should match the small-t form of ym(t) given by (|6.18fl . This 
requires 



(6.22) 



^As in 45 1, a < 6 (s -> 0+) 



as £ -> 0+ 
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Restoring the original variables and substituting w = u a results in a beta-function integral 
and the quantisation condition 



iB(4|)£w = (2fc + Z + §)7r, kZ. 



(6.23) 



Using (|4.16[) . the condition implied by (|6.2j) is recovered to next-to-leading order^. Notice 
that the full /-dependence has been recovered directly, without any need for resummation. 

Encouraged by this result, it is natural to see how much more of the asymptotic expansion 
found using the ODE/IM correspondence can be recovered by incorporating higher-order 
WKB corrections. Suppose that (|6.13l) has a solution of the form 



y(t) = exp 



1 oo 

-Tt n Sn(t) 



n=0 



(6.24) 



For (|6.13|) to be satisfied order-by-order in e, the derivatives S' n {t) must obey the following 
recursion relations: 



n-1 



S' (t) 



-VOW , 2S' S' n + £ S'jS'^ + S'U = (n > 1) 



(6.25) 



The first few terms of the solution are 
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(6.26) 



Keeping just 5o and S± constitutes the physical optics approximation employed above; further 
terms are very easily obtained using, for example, Mathematica. Near to zeros of Q(t) - the 
so-called turning points - the approximation breaks down and further work is needed to find 
the connection formulae for the continuation of WKB-like solutions of given order from one 
region of non- vanishing Q to another, just as was done above. Solutions found by continuation 
away from the two boundary conditions must then be matched to find the condition for Q to 
be such that the differential equation has an acceptable solution. This quickly becomes quite 
complicated, as exemplified by the calculations in section 10.7 of |l5] . 

In cases where Q is entire with just a pair of well-separated simple zeros on the real axis, 
Dunham [l^] found a remarkably simple formulation of the final condition, valid to all orders: 



1 r °° 

- i j^e n - 1 S' n {z)dz = 2kn, k£ 



(6.27) 



^Following Langer |42|. the same result can also be obtained starting from the approximation y(t) ~ 
l-Q(t)]- 1/4 C /2 (aJ-x(0 + l3Jx(0) in the region t w 0, with £ = \J ' —Q(u)du , and the boundary condition 
at t = requiring a — 0. 
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where the contour encloses the two turning points; it closes because for such a Q all of the 
functions S' n derived from (|6.25[) are either entire, or else have a pair of square root branch 
points which can be connected by a branch cut along the real axis. Notice that the contour 
can be taken to lie far from the two turning points, where the WKB series breaks down, 
and so there is no need to worry about connection formulae. All of the terms S' 2n+ i, n > 1, 
turn out to be total derivatives and can therefore be discarded, while 2 jS[ = —^Q'/Q and 
contributes a simple factor of ir/2 when integrated round the two zeros of Q, irrespective of 
any other details. Dunham's condition is therefore 



^e 2n - 1 S' 2n {z)dz = (2k+l)-n , keZ. (6.28) 

n=0 



However, this method is only directly relevant to the current problem if the angular mo- 
mentum is zero, and M is an integer. Extensive discussions of these cases can be found in |47j . 
and it is straightfowardly checked that the results found there match the expansion we ob- 
tained above using the ODE/IM correspondence. Note that for M £ Z there are no 'nonlocal' 
contributions to the asymptotic ()6.2j) . so the WKB series gives the complete answer. 

For more general cases there is only one simple turning point, the other being replaced 
by the singularity at z = 0, and the analysis just given does not apply. Nevertheless, the in- 
dependence of (|6.2j) together with the match with the results of |37j at special points suggests 
that the contributions to the asymptotic related to the local integrals of motion might still 
be obtained from a WKB series of the form (|6.28|) . suitably treated. The main difficulty is 
the fractional singularity at z = , which prevents the contour in Dunham's condition from 
closing. However, away from regions about z = and z = 1 which are vanishingly small as 
e — > 0, the all-orders WKB solution provides a good approximation to the true wavefunction. 
As an ad hoc measure, we replace the closed Dunham contour in (|6.27f) by a contour C which 
starts just below the origin, passes once round the turning point at z = 1, and returns to 
a point just above the origin. Since all of the terms S' 2n+1 , n > 1, are total derivatives of 
functions which are single- valued around z = 1, they again make no contribution even though 
the contour is no longer closed. Therefore we replace Dunham's condition (|6.28|) by 



JL 71=0 



s 2n ~ 1 S' 2n {z)dz = (2/fc + ^ + 2<5)vr , k£Z (6.29) 



where S' 2n is obtained from l)6.25j) using the Q given by ()6.14j) . and the factor 25 allows for 
a possible phase-shift caused by the singularity at the origin. From the initial asymptotic 
matching calculation, the leading (constant) part of 25 is A + 1/2, but further .E-dependent 
corrections can be expected - we shall return to this issue below. 

For concrete calculations, it is convenient to collapse the contour C onto real axis, and the 
square root singularities of the integrands at z = 1 then allow each integral to be replaced 
by 2 S' 2n (t) dt . A difficulty with this procedure is that the divergences in the WKB series 
at z = 1 are no longer avoided. We remedied this by multiplying each integrand S' 2n (t) by 
(i M / A - 1) K to force convergence, doing the definite integrals - still possible in closed form 
- and then setting k = at the end. It is a simple matter to mechanise the calculation 
with a few lines of Mathematica code, and, using (|6.4j) in reverse, the values of the local 
charges Ii, -Z3 and I§ given in (|5.24f) are easily reobtained; we also reproduced the formulae 
for I j and I9 quoted in |36| I39j , though these are rather too lengthy to be worth repeating 
here. It is interesting that the values of the conserved charges in the quantum field theory 
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can be recovered by such relatively-straightforward manipulations of the ordinary differential 
equation. 

Finally, the phase shift 5 should be analysed. We claim that this is related to the part 
of the asymptotic due to non-local conserved charges. To see why, focus on the behaviour 
near the origin by changing variables one more time in (|6.13j) . to u := e~ 2X t. The equation 
becomes 



d 2 



1 



du 2 4A 2 



1/A-2 



+ 



1 

4A 2 



£ 2M u (M+l)/X-2 



y{u) = . 



(6.30) 



Treating the final term as a perturbation results in a series of corrections to y{u) as powers 
of e 2M = E~^ I+l \ Comparing with l|fi.2j) . this is exactly the structure of the second series 
of terms emerging from the ODE/IM correspondence, encoding the values of the non-local 
conserved charges in the integrable model. 

In more detail, the perturbative treatment begins with the unperturbed equation 



d 2 



1 



.i/A-a 



y(u) = o. 



du 2 4A 2 
Two independent solutions j and n are 

j(u) = y/uJ\(u^) , Tl(u) = ■y/uYx(u^X) . 



(6.31) 



(6.32) 



Their Wronskian is W[j, n] := j{u)n' (u) — j' (u)n{u) = — v and using this fact it can be checked 
(as in for example |48l I49j ) that, with the boundary condition y(u) ~ j(u) as u — > 0, the full 
differential equation ()6.3U|) is equivalent to the integral equation 



y(u) =j(u) +e 



2M 



G(u\s) v(s)y(s) ds 



where 
and 



G(u\s) = vrA [n(u)j(s) - n(s)j(u)} 



v ( u ) = JL u (M+i)/A-2 
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(6.33) 



(6.34) 



(6.35) 



Formally, (|6.33|) can be solved by iteration. The zeroth order result is proportional to yjy(t), 
defined in 1)6.20(1 ; later terms give the Born series for the problem. 

To extract the corrected phase shift, the RHS of (|6.33|) can be rewritten as 



y(u)=j{u) l-e 2M 

The large-n asymptotics 

j{u) - 
then yield 

where 
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(6.38) 
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(6.39) 
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In problems with localised perturbing potentials, the integrals in . and also the function 

A(u) in (|6.38f) . tend to finite limits as u — > oo. Replacing y(s) on the RHS of l|6.39j) by the 
Born series generated by (j(i.33j) then gives the perturbative expansion of the phase shift. This 
must be treated with some caution here, since the perturbing potential l)6.35j) in general grows 
at infinity. However, at least at first order - the so-called Born approximation - it is possible 
to extract sensible results. 

The leading behaviour of (j6,33|) is simply y(u) ~ j(u). At this order there is no contribution 
from the denominator in (|6.39|) . and so the once-corrected phase shift is, formally, 

X 1 r°° 

= \ + \ + \e- {M+1) J™ x 2M+l Jx{xfdx + ... (6.40) 



The last integral converges for M + 1 + A > 0, 2M + 1 < , and in this region its value (found 
using Mathematica, or, for example, page 237 of [5U]) is 



M ; v / 2^r(-M)r(-M+A) v ; 

Using l|6.5|) . this reproduces the value of H\ given by Q5.25JI above. 

The agreement with previous results provides some retrospective justification for our pro- 
cedure, though we cannot rule out the appearance of subtleties at higher orders. Clearly, a 
more rigorous and systematic treatment would be desirable, especially given the difficulties in 
evaluating the eigenvalues of the nonlocal charges directly within the quantum field theory. 
However, we shall leave further investigation of this point for future work. 



7 Conclusions 

We have given a concrete application of the ODE/IM correspondence by showing how it can 
be used to obtain an analytic understanding of the level mergings in the model introduced 
by Bender and Boettcher, one of the longest-studied examples of VT symmetric quantum 
mechanics. The subtle mixture of perturbative and nonperturbative effects contributing to 
the recombination of levels in the generalisation of the model to include a centrifugal term 
is particularly striking, and shows once again the surprising richness of T>T symmetry as a 
source of interesting problems in mathematical physics. 

There are many questions left unanswered by our work, and we finish by mentioning just 
a few potentially interesting directions for further investigations. 

First of all, the origin of the perturbative and nonperturbative terms for M < 1 should be 
understood through more standard differential equation techniques. For the nonperturbative 
term responsible for the level-mergings, we expect that the complex WKB method will play a 
role, as in [S] . This could enable the study of models not treatable with the (less conventional) 
methods described in this paper. An alternative strategy for systems with potentials of the 
form P(x)/x 2 with P(x) a polynomial would be to use the set of functional relations (the 
fusion hierarchy) satisfied by the spectral determinant T(E), and this could also be explored. 

Our discussion of the nonperturbative term was restricted to M ~ 1; as mentioned above, 
when M goes beyond 1/2 or 2, the pole responsible for this term no longer crosses the steepest- 
descent contour and the asymptotic changes. This effect turns out to be most marked near 
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to M = 1/2 for I = 0, and in this region some further smoothing is needed before a good 
approximation for the energy levels can be obtained from the integral equation. Relevant 
techniques for analogous problems have been developed in |51| I52|. but we have yet to apply 
them to the current situation. It would be interesting to see whether the delicate asymptotic 
for the diverging ground state energy calculated in |^j could be recovered by such methods. 

The perturbative parts of the expansions for both the radial and the lateral problems 
are especially interesting for the ODE/IM correspondence, since they encode the values of 
conserved charges in the integrable quantum field theories. Our discussion of the relationship 
between the Born series and the nonlocal conserved charges at the end of section El was rather 
preliminary, and it would be nice to make the analysis fully rigorous and to push it to higher 
orders. In this respect, the sophisticated calculations of spectral zeta functions of Chudnovsky, 
Chudnovsky and Voros (SHI are likely to be relevant. The relationship between the local 
conserved charges and the WKB series is much clearer, but still needs to be understood on 
a more profound level. On the ODE side, efficient techniques have been developed for the 
evaluation of higher-order WKB contributions (see for example |54|). and there now seems to 
be scope to apply these methods to the study of integrable models. 

Finally, it would be worthwhile to extend these considerations to other integrable models, 
and other ordinary differential equations. Some results in this direction, for local conserved 
charges, can be found in [SHI ; two further examples to study would be the higher-order 
equations related to SU (n) Bethe ansatz systems discussed in [^ZJ , and the Schrodinger equa- 
tions for excited states found in |58| . 
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